Spin entanglement induced by spin-orbit interactions in coupled quantum dots 
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We theoretically explore the possibility of creating spin quantum entanglement in a system of two 
electrons confined respectively in two vertically coupled quantum dots in the presence of Rashba 
type spin-orbit coupling. We find that the system can be described by a generalized Jaynes - 
Cummings model of two modes bosons interacting with two spins. The lower excitation states of 
this model are calculated to reveal the underlying physics of the far infrared absorption spectra. The 
analytic perturbation approach shows that an effective transverse coupling of spins can be obtained 
by eliminating the orbital degrees of freedom in the large detuning limit. Here, the orbital degrees of 
freedom of the two electrons, which are described by two modes of bosons, serve as a quantized data 
bus to exchange the quantum information between two electrons. Then a nontrivial two-qubit logic 
gate is realized and spin entanglement between the two electrons is created by virtue of spin-orbit 
coupling. 

PACS numbers: 73.63.Kv, 03.65.-w, 03.67.Mn, 71.70.Ej 



I. INTRODUCTION 

Control and manipulation of the spin degree of free- 
dom become one of the most important topics both in 
spintronicsl'l and in quantum information processing 0. 
The spin-orbit interactions (SOI) in semiconductor het- 
erostructures provides the ways to couple spin with spa- 
tial degree of freedom, and consequently has attracted 
more and more attention in recent years. The spin prop- 
erties of a few electrons confined in semiconductor nanos- 
tructures, such as quantum dots^, H, |M IE S cou- 
pled quantum dots^], quantum rings [l(l| and quantum 
wires have been studied. The results show that the 
carrier spin properties are strongly affected in the pres- 
ence of the SOI, and novel features emerge in these nonos- 
tructures compared with the traditional ones without 
SOI. 

On the other hand, spin confined in quantum dot is a 
natural choice for the physical realization of qubit. This 
kind of system is considered as an important candidate 
for solid state based quantum computing. Among vari- 
ous approaches to implement quantum information pro- 
cessing using quantum dot systems, the optical method, 
including the classical laser field[l3 and the quantized 
cavity modes [HI, has been proposed to create entan- 
glement and to realize the single and double qubit logic 
gate. In Ref. 13, a scheme with quantum dots embed- 
ded in an optical cavity was designed, so that the cavity 
mode can serve as a data bus and induce a spin-spin in- 
teraction. This kind of cavity-mediated two-qubit gate 
is studied in details for several other solid state systems 
very recentlvfl^. 

Now we notice that the SOI phenomena in nanostruc- 
tures were investigated with the help of the quantum 
optics method. Taking advantage of the tunability of the 
SOI strength, an experiment to observe coherent oscilla- 
tions in a single quantum dot was proposed in ref. 3 • 
this proposal, the orbital degrees of freedom are modeled 
by two boson modes. Under the rotating wave approx- 



imation (RWA), the SOI of the electron confined in the 
single quantum dot, is reduced to a Jaynes- Cummings 
(JC) model, which is very typical model in quantum op- 
tics. This analogy between the SOI in a quantum dot and 
the JC model in cavity QED suggests us that it is possi- 
ble to make use of the orbital degree of freedom, instead 
of the real optical cavity modes, to serve as a quantized 
data bus and then to induce a spin entanglement . 

In this paper, we propose and study a model of two 
electrons confined in two quantum dots respectively with 
Rashba type SOI, and explore the possibility of realiz- 
ing a two-qubit logic gate or creating spin entanglement 
with this system. For simplicity, we consider two ver- 
tically coupled quantum dots (CQDs) with two dimen- 
sional parabolic confinements. In the case with strong 
confinement and large interdot separation, the Coulomb 
interaction between the two electrons is approximately 
expanded in a quadratic form, and then the orbital part 
of the two electron system can be reduced to four-mode 
bosons. Under the RWA, only two of the four modes 
are coupled with the spin degrees of freedom. Then it is 
proved that the total system with SOI can be effectively 
described by a generalized JC model with the coupling 
between two bosons and two spins. By diagonalizing the 
Hamiltonian directly in the lower excitation subspace, 
the eigenvalues and the corresponding eigenfunction are 
obtained exactly. The far infrared (FIR) absorption spec- 
tra are calculated according to these analytical solutions, 
which help us to understand the underlying physics of the 
spectra. To get an effective Hamiltonian of spin-spin in- 
teraction, we perform the Frohlich transformation in the 
large detuning limit. This effective Hamiltonian can dy- 
namically drive a two-qubit logic gate operation. By us- 
ing the conventional material parameters, our numerical 
estimation shows that the effective spin interaction in- 
duced by SOI is strong enough, compared to the spin de- 
coherence in low dimensional semiconductor structures. 
It is feasible experimentally to implement a two-qubit 
logic gate and thus produce quantum entanglement. 
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FIG. 1: Schematic illustration of the vertically stacked quan- 
tum dots. Each quantum dot confined an electron. We as- 
sume that the interdot separation zq is larger than the char- 
acteristic length of quantum dot confinement la- 



the total Hamiltonian reads: 
^ = E f ^ + ^ (^^) + ^sl + H^z^) + ycoui. (4) 



Here, — e , /is and Sq are the electron charge, Bohr mag- 
neton and dielectric constant in vacuum, mo, g and e are 
the material related parameters of effective mass, Lan- 
dre g- factor and the relative dielectric constant, respec- 
tively. Tli = Pi + eA(ri) is the canonical momentum 
and A(ri) = B {—yi/2, Xi/2, 0) is the vector potential for 
magnetic field B — Bz. 

In order to simplify the Coulomb interaction, we con- 
sider a special case in which the interdot separation zq 
is much larger than the lateral confinement characteris- 
tic length Iq = y^h/moU!o, i.e. {Iq/zq) << 1. Then 
we expand the Coulomb interaction as a power series of 
the relative coordinate r = |ri — up to the second 
order 113 [3: 

Vcoui (r) ^ Vo - ^moujfr^. (5) 



The paper is organized as follows. In Section II, the 
two modes JC model is derived from the conventional 
Hamiltonian of two electrons confined in two quantum 
dots. Analytical solution of lower excitation states and 
FIR spectra are presented in Section III. In Section IV, 
we show that the perturbation treatment gives the effec- 
tive transverse spin-spin interaction in the large detuning 
limit, and demonstrate that a two-qubit logic gate and 
quantum entanglement can be achieved in this kind sys- 
tem with SOI. 



II. TWO MODES JC MODEL FOR 
VERTICALLY CQDS WITH SOI 

We consider two vertically CQDs with an interdot sep- 
aration Zq, see Fig^ Each quantum dot is described 
by a two dimensional parabolic confinement potential 
V (vi) = moujQrf/2 with the basic frequency ujq. Includ- 



ing the Rashba type SOI 
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the Zeeman term 
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and the Coulomb interaction 
Kouzdri - r2|) = 



47re£o ri - r2 



(1) 



(2) 



(3) 



Here, we have defined Vq — e^/ [A.-keeqZq), and hbj\ = 



We also assume that electrons are strictly 



confined in each quantum dots planes, and then neglect 
the overlap of their wavefunctions. 

In the center of mass (CM) reference of frame defined 
by R = (ri -|- /2 and r = ri — r2,we have the CM and 
relative momentums, and the corresponding angular mo- 
mentums P = pi -I- p2, p = ^ (pi — P2) , L = R x P, and 
1 = r X p, where M — 2mo,m — mo/2,r is the relative 
coordinate and R the CM coordinate. The orbital part 
of the Hamiltonian is expressed in a quadrature form of 
these coordinates: 



E l^ + ^i^^) ]+Veoul 



\2too 



2M^2 

1 2 2 1 , 

1 — muj r H — oJAz 

2m 2 2 



(6) 



Here, the cyclone frequency is lOc = eB/mQ, and fre- 
quency of CM and relative motion fl 
respectively. 



+ and 



u! — •\/rP^^2a^, respectively. Note that the effect of 
the Coulomb repulsion is reducing the relative motion to 
a lower frequency compared to the CM motion. In our 
model the requirement {Iq/zq)^ << 1 ensures that the 
w > is satisfied even when B ~ 0. 

We define the boson operators a^:, a+, ax, at- of x com- 
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Quantity Value 


Unit 


Quantity Value Unit 


mo 0.042 






20 meV 


e 14.6 






9.5 nm 


9 -14 




zo 


20 nm 


Q 10 


meV-nm 







TABLE I: Parameters used in the calculations 



ponents by 




and the boson operators ay,a^ ^ay, fly of y components 
are defined in the same way. Let 





{ax ^ 


' iay) /V2, 




a — 


{ax + 


iay) /V2, 




B = 


{ax - 


- iay) /V2, 


(8) 




{ax - 


iay) /V2. 





Then we rewrite orbit Hamiltonian ^ 

i/orbit = {A+A + 1/2) + hujB {B+B + 1/2) 

+ huja {a+a + 1/2) + hcji, {b+b + 1/2) (9) 

in terms of four mode bosons A, a, i3,and b, with their 
frequencies respectively 

OJa^^ - itJc, UJb ^UJc, (10) 

1 1 

= - -Wc, Wft = + -Wc (11) 

These four frequencies, together with the Zeeman energy 
huJz = \g\^BB are drawn in Fig|21 with respect to the 
magnetic field B. Note that the Zeeman energy huJz 
and the boson frequencies uja and Ua reach the resonant 
regime at i? ~ IIT with the parameters listed in Table 
n We also draw the lower energy spectra of the orbit 
Hamiltonian ^ in Figl^Jb). We will focus on how these 
states are affected in the presence of SOI in the following 
sections. 

In terms of the four boson operators defined by Eqs. 
(0 and JHJ, the SOI Hamiltonian Hso, after some 
straightforward algebra, can be rewritten as: 

HiS'^ = gA-A{a,+ +<j,+ ) 

+ 9a- a {ai+ - (72+) + h.c. (12) 




'o 3 6 9 12 15 

Magnetic Field B (T) 



FIG. 2: (a) Energy dispersions with respect to the magnetic 
field B of the four modes of the boson frequencies uja (red 
solid line), uJa (blue solid line), ujb (red dashed line), and uji, 
(blue dashed line), and the Zeeman energy Wz (black line), 
(b) Energy spectrum of Hq. Different lines correspond to 
diff'erent spin and orbit states, which are indicated explicitly 
in the figure. The parameters used in calculation are listed in 
Table El 

Note that, due to the negative value of the Landre g- 
factor, a unitary rotation CTz i— > —cTz and cr± i— > —crip has 
been performed. To obtain the interaction Hamiltonian 
above, we have used the RWA to neglect the counter- 
rotating terms like a+B~^ , <^+b~^ , cr_i?,and cr_6. This ap- 
proximation has been verified numerically for single elec- 
tron case in ReL The coupling strengths qa and ga 
are defined as follows: 

-"V¥('-S) <-) 

So far, we have obtained a generalized JC model where 
two- mode bosons interact with two spins. In the follow- 
ing sections we will further demonstrate how the orbit 
motion induce a spin-spin entanglement in the presence 
of SOI. 



III. LOWER EXCITATION STATES AND FIR 
SPECTRA 

In this section, we calculate the eigenenergies and 
eigenstates of the low excitation states. Notice that the 
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FIG. 3: (a) Exact energy spectra given by Egs. 11231 and 11241 . 
The arrows indicate the dipole allowed transition from \GS} 
(black line) to 1$^) (red lines). |$J) (blue lines) is dipole 
inactive, (b) FIR spectra for increasing magnetic field B. 
The anti-crossing is clearly shown. A Lorentz profile function 
T/tv[{u> — u)fi) — T^] is used to replace the 5 function. The 
phenomenological broaden factor F = 0.2meV. Spectra are 
normalized to their maxima and offset for clarity. 



this subspace. We define 
1 



\To) 



V2 
1 



(|OA,Oa,T,i)-|OA,Oa,i,T)), 
(|OA,Oa,T,i) + |OA,Oa,i,T)), (18) 



|T^l) = |U,Oa,i,i),|rfi) = |OA,la,ia). 

The eigenstates of the Hamiltonian in this subspace are 





= sin 


Oa 
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+ COS 




(19) 




= cos 


Oa 
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— sin 




(20) 




— sin 






-f cos 


\\S). 


(21) 




= cos 


Oa 

2 




— sin 




(22) 



and the corresponding eigenvalues are: 




(23) 
(24) 



Here, the zero point energies of the four boson modes are 
omitted. Oam is defined as 



tan( 



(25) 



total excitation number operator 



N ^a+a + A+A+-{ai,+a2,) (15) 



commutes with H = i/orbit + H^z + + H^o"^- For 
a given integer N , which is the eigenvalue of iV, the di- 
mension of the invariant subspace y(^) is 4iV + 4 for 
TV > 0. 

The lowest subspace V^~^\ corresponding to iV = — 1, 
is of one dimension. The the ground state can be directly 
written as \GS) = |0a, Oa, Ob, Ob, |, |), which means the 
excitation numbers of boson modes 



where 9A,a & [0,7r], and Aam =^z- ^A,a- 

With the help of these exact eigenenergies and the 
eigenstates obtained above, we can calculate the FIR ab- 
sorption spectra analytically. To this end, we consider 
the time-dependent Hamiltonian of the system in an op- 
tical field is 



il{t) + H'e' 



(26) 



where H' is the time dependent term induced by the 
classical optical field 



^ — ' L m 

1=1.2 



aea , 

— (a, X e)^ (27) 



A+A = a+a = B+B = 6+6 = 



(16) 



and both spins are in the down state. The corresponding 
eigenenergy is 



Egs = 2 i^A +i^B + (^a + UJb- ^^z) ■ (17) 



The second lowest subspace which is of four di- 

mension. The Hamiltonian can be exactly diagnolized in 



where a is the radiation field amplitude, and e ~ 
{x — iy) I \/2 is the polarization vectors for circular po- 
larized light. 

The absorption coefficient is calculated according to 
the Fermi golden rule[l9t[20l|: 



a(^)«c.^|(/|ff'|z)|'5(o 
/ 



(28) 



Here, Sw/i = Ej — Ei, \i) = \GS) is the ground state, 
and I/) stands for the excited states. We will focus on 
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the four lowest excited states and \^^)- According 

to Eqs. {TI) and the perturbation H' is obtained as 

(29) 

Note that we have omitted the terms related to the 
B modes, which is not involved in the initial and final 
states we are considering and thus does not contribute 
to the absorption coefficient. Our analytical results give 
the FIR spectra of obvious physical meanings. From 
Ea. (|29|l above, we find that the first term is spin in- 
dependent, and it provides an CM angular excitation. 
This CM angular excitation, which is the consequence of 
Kohn theoremj^J, exists even in the absence of SOI. The 
second term, which is spin dependent, is due to the pres- 
ence of SOI. This term contributes an excitation of the 
two spins symmetrically, i.e. ||J,) i — > |ti) + lit)- Thus, 
H' only couples the ground state \GS) with the CM ex- 
citation states 1$^), and are inactive in this case. 
The matrix elements of {^^\H'\GS) are calculated, and 
the FIR spectra are shown in FigOfb) . 



IV. QUANTUM ENTANGLEMENT IN LARGE 
DETUNING LIMIT 

Due to the linearly increasing of the dimension of the 
invariance subspace y(^) with respect to N, it is difficult 



to obtain a compact solution of the eigenvalue problem 
for > 1. Instead of the exact solution, the approx- 
imate solution by perturbation theory will be given in 
this section. In the perturbation available regime, we de- 
rive an effective transverse spin-spin interaction Hamil- 
tonian. This Hamiltonian can induce a two-qubit logic 
gate, and can be used to produce a controllable quantum 
entanglement. 

We summarize the Hamiltonian obtained 

Ho = Horhit + ^fi^z {(^iz +cr2z) , (30) 

Hi = gA ■ A{ai+ + (J2+) 

+ ga ■ a {ai+ - a2+) + h.c. (31) 

and then consider its reduction in the large detuning 
limit, i.e. AA,a >> gA,a- In this limit, we perform the 
Frohlich transform with the operator 

S = (ai_ + <72-) + (ai_ - f72-)) - h.c. 

(32) 

and the effective Hamiltonian exp(— 5*)^^ exp(S') is calcu- 
lated up to the second order as: 

Hs^Ho + ^[HuS]. (33) 

Here, the second term in the r.h.s can be written explic- 
itly 



- [Hi,S] = h4 (cri+cr2- + cri_cr2+) 



9Aga 



<y2z 



1^ {A+A + 1/2) + 1^ {a+a + 1/2) 
Aa Aa 



{<Jlz 



9\ 



gl 



(34) 



Note that the first term of Ea. l34|l is the effective trans- 
verse spin-spin coupling induced by SOI. The effective 
coupling strength 



eigenenergies 



n£, = g\IAA - gl/Aa 



(35) 



depends on (i) the SOI strength a (see Eas. ljTHl and 
(llSfl ). and (ii) the frequency difference between loa and 
LOa, which is the consequence of the Coulomb interaction. 
We note that, in the subspace V^'^\ the second term of 
Ea. (|34|l vanishes, and the remaining terms commute with 
the total spin 0-2 = (o-i -ho-a) . Thus we can denote 
the eigenstates by spin singlet \S) = (|ti) — lit)) /\/2 
and triplet |To) = (|ti) + UT)) /\/2, |Ti) = |tT), and 
|T_i) = III) in the large detuning limit in subspace 
V'^^\ Diagonalizing the Hamiltonian (|34|l . we obtain the 







25i 
Aa 


^|T»,) 


= -A,- 


-^gl 

Aa' 




A^' 




E\s) 


A, 





(36) 
(37) 
(38) 
(39) 



for the different spin states corresponding to the exact 
solutions in Ea. (|24|l . Here, we also omit the zero point 
energies. 

The above energies are drawn as the function of mag- 
netic field in Fig^c) in comparison with the exact solu- 
tion. The exact consistency of the two solutions in the 
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FIG. 4: (a) The ratio qa/^a (black line) and, Qa/^a (red 
line). The perturbation treatment is valid when both qa/ ^a 
and gal ^a. less than unity, (b) The effective transverse 
spin coupling strength induced by SOC. (c) A comparison 
of eigenenergies between exact (Eos. 1)23^ and (|2U, lines) and 
perturbative(Eqs.l(3ni ~ H-ffl|l ■ scattered dots) solutions in the 
large detuning regime. Our perturbation treatment is valid 
in the unshaded region. 



FIG. 5: (a) Probability Pi (t) as a function of t. Black, red, 
blue, and orange solid lines represent the probabilities of 
\n,OA, Oa), UT,OA,Oa), Ui,lA,Oa), aud Ui,0.4,la) respec- 
tively. The dashed lines are the probabilities generated by 
the standard V zSWAP gate. (b)Gate fidelity as a function of 
t. The black line indicates the ideal V iSWAP operation. The 
blue and red lines represent the fidelity function generated by 
the second order approximation Hamiltonian and the original 
Hamiltonian (see text) respectively. 



large detuning regime confirms the validity of our per- 
turbation treatment. Fig^a) shows effective spin cou- 
pling strength induced by SOI in our model. In the per- 
turbation valid regime, for example a,t B = lOT (the 
point denoted by X), we have the spin coupling strength 
~ 20/ieV^. This value is comparable with that in the 
proposal of the spin-spin coupling induced by the by elec- 
tromagnetic field in cavitvflSj. Thus, it is feasible to 
realize the two-qubit gate operation during the long co- 
herence time of conduction band electrons. 

The spin transverse coupling described by t he first 
term in the r.h.s. of Ea. (|M|l generates an ideal V iSWAP 
gate|220| at a specific time to = 7r/4^. On the other 
hand, it is worthy to notice that the unitary transfor- 
mation exp(— 5")^? exp(S') may induce an excitation of 
the boson modes, and then causes gate error during the 
time evolution of the spin states. To explicitly show the 
unwanted boson modes excitation and to examine the re- 
liability of our model, we consider the time evolution of 
the initial state (0)) = |iT,OA,Oa): 



|V(t)) =t/(t)|^(0)). 



where 



Notice that |-0 (0)) belongs to the subspace V'-^^ and the 
time evolution will be restricted in this subspace. Fi^^Ia) 
gives the probabilities Pi {t) = | (t)) | of finding state 
\i) at time t, where |z) stands for the four bases of the 
subspace 1/(0): |n,OA,Oa), \ltOA,Oa), Ui,U,Oa), and 
|ii,OA,lo) respectively. Besides the states |ti,OA,Oa) 
and lltjOyijOa), the unwanted boson modes excitation 
states III, 1a, Oq), and \ll, Oa, la) also have nonzero pop- 
ulations. This populations will induce the gate error. 
Furthermore, we calculate the fidelity function defined 
by 



Fit)^ V(0) C/, 



%/iSWAP 



[/(<) 7^(0) 



(40) 



U (t) = exp [it (e'^Hse^'^) /h] = e"^ exp (itHs/h) e~ 



where U ^ .^swap ideal ViSWAP gate operator. We 

notice that the fidelity function F (t) reaches its maxi- 
mum slightly less than unity at t = to, and the high- 
frequency oscillations appear due to the boson modes ex- 
citation mentioned above. We also examine the fidelity 
function generated by the original Hamiltonian (|31|l . i.e. 
U (t) = exp [—it {Hq + Hi) /H]. In this case the original 
Hamiltonian H31|l gives a high fidelity F {t') ~ 1 at a 
different time t' ~ 1.2to, which can be regarded as the 
higher order correction in comparison to the approximate 
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Hamiltonian (|33|l . Finally, Fig|Sl shows that spin entan- 
glement can be created in this system by adjusting the 
operation time. 

V. CONCLUSION 

In this paper, we have considered a system of two ver- 
tically CQDs each containing an electron in the presence 
of Rashba type SOI. We theoretically demonstrate that 
it is possible to create spin entanglement in this kind of 
system by using the SO coupling. With the large inter- 
dot separation case, the Coulomb interaction between the 
two electrons is approximately expressed in a quadratic 
form. And then two-boson-two-spin interacting model is 
derived in the RWA from the original Hamiltonian. We 
give the exact solution of the low excitation states ana- 
lytically. This solution helps us reveal the physics under 
the FIR spectra near the resonant point. Perturbation 
treatment in the large detuning case shows that, similar 
to the quantum dot embedded in an optical cavity, the 
orbital freedoms play a role of quantized data bus via the 



Coulomb interaction and SOI in this system. An effec- 
tive Hamiltonian of spin-spin interaction is obtained in 
the perturbation regime by eliminating the orbital free- 
doms. This Hamiltonian provides a tow-qubit operation, 
which is essential in quantum information processing. 

Finally, we would like to point out that using the ef- 
fective inter-spin coupling to create spin entanglement is 
feasible to be controlled and measured. From the dis- 
cussion above, we know that the effective spin coupling 
strength can be controlled by external magnetic field. On 
the other hand, the tuneable strength of SOI a 
principle, also enable us to switch on and off the effective 
inter-spin coupling by external gates conveniently. To 
probe the quantum entanglement of spin system, simi- 
lar method in the protocol proposed in ref. can be 
used, where the information stored in the spin degrees of 
freedom is converted to the charge states, and then the 
charge states can be detected. 

This work is funded by NSFC with grant Nos. 
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NFRPC with Nos. 2001CB309310, 2005CB724508, 
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